Quantum and classical thermal correlations in the XY spin-1/2 chain 



J. Maziero,^ H. C. Guzman,^ L. C. Celeri/ M. S. Sarandy,^ and R. M. Serra^ 

^Centra de Ciencias Naturais e Humanas, Universidade Federal do ABC, 
R. Santa Adelia 166, 09210-170, Santo Andre, Sao Paulo, Brazil 
^Instituto de Fisica, Universidade Federal Fluminense, 
Av. Gal. Milton Tavares de Souza s/n, Gragoatd, 24210-346, Niteroi, Rio de Janeiro, Brazil. 

We investigate pairwise quantum correlation as measured by the quantum discord as well as its 
classical counterpart in the thermodynamic limit of anisotropic XY spin-1/2 chains in a transverse 
magnetic field for both zero and finite temperatures. Analytical expressions for both classical and 
quantum correlations are obtained for spin pairs at any distance. In the case of zero temperature, it 
is shown that the quantum discord for spin pairs farther than second-neighbors is able to characterize 
a quantum phase transition, even though pairwise entanglement is absent for such distances. For 
finite temperatures, we show that quantum correlations can be increased with temperature in the 
presence of a magnetic field. Moreover, in the XX limit, the thermal quantum discord is found to 
be dominant over classical correlation while the opposite scenario takes place for the transverse field 
Ising model limit. 
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I. INTRODUCTION 

Since its conception, quantum mechanics has been ex- 
tensively applied to condensed matter systems [l|. This 
endeavour resulted in important conceptual and tech- 
nological advances. In the last few decades, this sce- 
nario has gained even more strength with the born of 
quantum information science (QIS) [2]. In this context, 
quantum spin systems have played a central role in sev- 
eral applications for both quantum communication [3] 
and quantum computation [4]. Indeed, quantum spins 
models describe the effective interactions in a variety of 
physical systems [5], e.g., quantum Hall systems, high- 
temperature superconductors, heavy fermions, and mag- 
netic compounds. More generally, low dimensional sys- 
tems as those described by interacting spins are partic- 
ularly interesting because of the presence of typically 
pronounced quantum fluctuations as well as the possi- 
bility of their realization by several distinct physical ap- 
proaches [l, 0- 

A key concept in QIS is the quantum correlation 
among parts of a composite quantum system, which is 
a fundamental resource for several applications in quan- 
tum information [2j. The existence of quantum corre- 
lations were first noted in nonseparable (i.e., entangled) 
states. Entanglement was pointed out by Schrodinger in 
1935 [H as the characteristic trait of quantum mechan- 
ics. Since then, entanglement was basically the only kind 
of quantum correlation theoretically and experimentally 
explored. However, in the past few years, it has been 
realized that there exist nonclassical correlations which 
are not captured by entanglement measures [9*, To]. In 
general, a quantum correlation, which can be measured 
by the quantum discord (QD) [9], often arises as a con- 
sequence of coherence between different partitions in a 
quantum system, being present even in separable states. 
Recently, QD was analyzed in a number of contexts, e.g., 
low dimensional spin models pj]-[l5| , open quantum sys- 



tems p!6l-[T9^, biological [20], and relativistic [21^ systems. 
Moreover, there exist strong indications that the QD is 
the resource responsible for the speed up in the model of 
computation known as deterministic quantum computa- 
tion with one quantum bit j22|, [23| . 

In this article we consider pairwise QD in an infinite 
anisotropic XY spin-1/2 chain in the presence of an exter- 
nal transverse magnetic field for both zero and finite tem- 
peratures. Our aim is to explore, in the thermodynami- 
cal regime, the behavior of QD for spin pairs arbitrarily 
distant and also take into account the effect of tempera- 
ture on the behavior of correlations. Such contributions, 
which have not been considered in previous works, will 
be shown to bring several new effects to the subject. In 
particular, as we will show, the QD for spin pairs more 
distant than second-neighbors is able to characterize a 
quantum phase transition (QPT). This is a remarkable 
behavior, since a signature of the QPT is then available 
even for distances where pairwise entanglement is ab- 
sent. Moreover, we will show that the QD may increase 
with both temperature and magnetic field for certain re- 
gions of parameter space. This result will extend, to the 
thermodynamic limit, the previous analysis for two-spin 
Hamiltonians reported in Ref. [l3[. Finally, we will dis- 
cuss the dominance of quantum correlation over classical 
correlation for different limits of the XY model, showing 
that the QD is greater than its classical counterpart for 
the isotropic limit (XX model), with the opposite sce- 
nario taking place for the transverse field Ising model. 
These results generalize those of Refs. [ll| and |12| . 



II. CORRELATIONS IN BIPARTITE 
QUANTUM SYSTEMS 

The information-theoretical measure of the total cor- 
relation between the partitions of a bipartite quantum 
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state pAB is the quantum mutual information [2J, |25| 

I{pab) = S{pa) + S{pb)-S{pab), (1) 

where S {p) = — Trplog2P is the von Neumann entropy 
and Pa {pb) is the reduced- density operator of the parti- 
tion A {B). The nonclassical correlation present in pab 
can be quantified by the quantum discord, which is de- 
fined as [9] 



D (pab) = I (pab) - C (pab) , 



(2) 



where 



C (pab) = S (pa) - min S^u,} {pa\b ) (3) 



is the classical correlation present in the composite 
state Pab [26]. In Eq. (j3j), the conditional entropy 
'^{Uj} {pa\b) can be defined as 



^{n,} {pa\b) =^QjS{p\), 



(4) 



with Qj = Tr [(1a (8) n^ ) Pab (1a (8) n^ )] and p\ = 
TTB[{lA'S)Ilj) Pab {'^A'S)Ilj)]/ Qj. The minimum in 
Eq. ([3]) is taken over a complete set of projective mea- 
sures {Hj} on the partition B. For pure states, we have 
that both quantum and classical correlations are equal 
to entanglement entropy [l^, [HI- On the other hand, for 
mixed states, entanglement is only a part of this nonclas- 
sical correlation |i,l2i,[23j. 

In order to compare our results for the QD with those 
for pairwise entanglement, we will use the entanglement 
of formation (EoF) as a measure of entanglement. The 
concurrence (cone) is monotonically related to the entan- 
glement of formation by the following expression 



EoF{pab) = y^bi 



1 



^1 - (conc(pAB))^ y/ 2| , 



where 1-Lbin{x) = — xlog2X — (1 — x)log2(l — x) is the 
binary entropy. For two qubits the concurrence reduces 
to the simple form [27] 



cone (pab) = max ^0, ^/Xl 



(5) 

where {i = 1,2,3,4) are the eigenvalues of PabPab 

in decreasing order, pAB = {c^a ^ ^b) Pab (^a ^ 
with p\^ being the conjugate of pab in any basis, and 
cr^ is the ^/-component of the spin- 1/2 Pauli operator 
for the partition a {a = A^B). We note that all the 
correlations quantifiers aforementioned are measured in 
bits, as usual. 



III. THERMAL CORRELATIONS IN THE 
ANISOTROPIC XY SPIN CHAIN 

The ID XY model in a transverse field describes a 
chain of spins anisotropically interacting in the xy spin 



plane under the effect of a magnetic field in the z direc- 
tion. The system is governed by the Hamiltonian 



H 



N-l f 

j=0 



[(l + 7)^J^J+i + (l-7)^^. 



y^y 1 



(6) 



where {k = x^y, z) is the k component of the spin- 
1/2 Pauli operator acting on site j state space, 7 is the 
degree of anisotropy (we take for simplicity < 7 < 1), 
and A provides the strength of the inverse of the external 
transverse magnetic field. We will be interested in the 
limit of an infinite chain, namely, N ^ 00. 

The XY model is exactly solvable [H, [29j. The 
Hamiltonian can be diagonalized via a Jordan- Wigner 
map followed by a Bogoliubov transformation (see, e.g., 
Ref. [30]). By taking the thermal ground state, the re- 
duced density operator for the sites and n reads [3l| 



P07 



= \ I Ion + {ct') (ag + O + ^(a^a^)a^a^ 

I k=l 



(7) 

with /on being the identity operator acting on the state 
space of the sites On. Although an unbroken state, pon is 
able to provide an exact description of the critical behav- 
ior as well as its scaling in finite systems [3ll-[35i (For a 
detailed treatment of the spontaneous symmetry break- 
ing at zero temperature see Refs. [36"-^38]). 

Since the system is invariant by translations, the ele- 
ments of the two site reduced-density operator depends 
only on the distance (n) between the sites. The trans- 
verse magnetization is given by [28] 

+ A cos (j)) tanh (/^a;^) 



Jo 



-d^. (8) 



where cj^ = (7Asin(/))^ + (1 + Acos(/))^/2 and (3 = 
1/kT with k being the Boltzmann's constant and T the 
absolute temperature. The two-point correlation func- 
tions read [29j 
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The total correlation in ^ is quantified by the 
quantum mutual information as I{pon) = S{po) + 

S{pn) - S{pon) with S{po) = S{pn) = -J2i=o{l'^ + 

(-l)^^(cT^)]/2}log2{[l + (-l)*(a-)]/2} and 5(pon) = 
-El,j= o(gi log2 + 6 log2 Ci), wher e = [1 + (agO + 
(-1)V(K<T-) - (<ToV?())2 + 4(<T-)2]/4 and = [1 - 
('^oO + (-l)^({<T^O + {<t^<tJ1))]/4. We can com- 
pute the QD and its classical counterpart by extremiz- 
ing Eqs. (|2]) and (|3]) over the following complete set of 
orthonormal projectors {11^ = \Q(3) (O/3I , /3 =||,±} onto 
the nth nearest-neighbor, where |0||) = cos(^/2) |0)^ + 
e^^sin(l9/2)|l)^ and IGx) = e"^^ sin((9/2) |0)^ - 
cos(^/2) Remarkably, a numerical analysis im- 

plies that the extremization is achieved, for any values 
of 7, A, and T, by the choice = 7r/2 and cp = 
[39]. Then, the relevant measurement for the model 
is given by {|+)(+U-)(-|}, with |+) and |-) denot- 
ing the up and down spins in the x direction, namely, 
|±) = (I t)"^ ± I |)^)/V^. This result generalizes the ex- 
tremization obtained for the transverse field Ising model 
at zero temperature in Ref. [12], and allows us to write 
the classical correlation as 

C(pon) = y^bin (Pl) - y^bin {P2) , (13) 

where 

pi = i (1 + (<T^)) , (14a) 




(14b) 



Thus the quantum correlation in state (O is simply given 

by 

D{pon)= I{pOn)- C(pon). (15) 

A. Correlations at zero temperature and QPTs 

Let us first consider the XY model at zero temperature. 
Such a model has a quantum phase diagram displayed in 
Fig. 1 (see, e.g., [28, 29, 41]). 

The QD for the first, second, third, and fourth nearest- 
neighbors in the thermal ground state ([7]), close to zero 
temperature, is displayed in Fig. 2 (a)-(d). As ex- 
pected, QD decreases as we increase the distance be- 
tween the sites. However, we also see a clear difference 
in the amount of quantum correlation between the re- 
gions where A < 1 and A > 1. Although the maximum 
value of QD decreases as we increase the distance be- 
tween the sites, the slope in the critical region gets more 
evident for far neighbors. The maximum increasing rate 
of the QD as function of A occurs at the quantum phase 
transition line (A = 1). We also note that nonclassical 
correlation is created when the magnetic field increases. 
The derivative of the quantum discord between fourth 
nearest-neighbors with respect to A is depicted in Fig. 3. 
So, the quantum discord between far neighbors can be 
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FIG. 1: (Color online) Quantum phase diagram for the 
anisotropic XY spin 1/2 chain. The XX model obtained by 
setting 7 = 0, displays a critical line for A G [0, 1]. The Ising 
model obtained for 7 = 1, exhibits a critical point at A = 1. 



used to characterize the QPT. It is important to mention 
that, even when entanglement is not present, the QPT 
can be clearly revealed through the singular behavior of 
QD. In Fig. 2 (e) and (f) we note that, for a consider- 
able range of 7, EoF is zero (there is no entanglement) 
for third and fourth nearest-neighbors. In the particular 
case of the transverse field Ising model (7 = 1), entangle- 
ment is indeed completely absent for sites farther than 
second nearest-neighbors [32|- On the other hand, QD is 
non-null (where the entanglement is zero) and its behav- 
ior reveals that, for 7 = 1 and A > 1, we have long-range 
nonclassical correlation between far sites. It is worth- 
while to observe that the pairwise classical correlation 
can also be employed to detect a QPT in a very evident 
way. In the Ising model, for example, the first derivative 
of the classical correlation with respect to A is not ana- 
lytic at the critical point, while the second derivative of 
quantum discord presents such a nonanalicity [12]. This 
behavior of the classical correlation also holds in the XY 
model for the whole range of values of 7 considered here, 
as depicted in Fig. 4. 



B. Correlations at finite temperatures 

In order to introduce finite temperature, we begin by 
considering the thermal state of the XX model (7 = 0). 
For this limit, we plot the thermal classical and quan- 
tum correlations for the second nearest-neighbors in Fig. 
5. Similar results are obtained for others than second 
nearest-neighbors. We note that the quantum discord is 
typically greater than its classical counterpart and it in- 
deed may increase with temperature (in a given region) 
for some values of A. It is evidenced near to the critical 
value A = 1. This result extends, to the thermodynamic 
limit, the previous analysis for two-spins Hamiltonians 
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FIG. 2: Quantum discord between (a) first, (b) second, (c) 
third, and (d) fourth nearest-neighbors and EoF between 
(e) third and (f) fourth nearest-neighbors as a function of 
anisotropy (7) and A at zero temperature. 
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FIG. 3: (Color online) Derivative of the quantum discord be- 
tween fourth nearest-neighbors with respect to A at zero tem- 
perature. 



FIG. 4: (Color online) Derivative of the classical correlation 
between fourth nearest-neighbors with respect to A at zero 
temperature. 



reported in Ref. [l3| . The increasing of QD with temper- 
ature when a magnetic field is present is a consequence of 
the fact that, when the field is turned on, the ground state 
tends to be less correlated than some low-lying excited 
states. So the effect of the temperature is to populate 
such correlated excited states, leading to the net effect of 
an increasing of QD. Naturally, this effect tends to dis- 
appear as the temperature gets too large. This is similar 
to the behavior of entanglement observed in Ref. [42j . 
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FIG. 5: (a) Classical and (b) quantum correlations between 
second nearest-neighbors in the XX model as a function of 
temperature (kT) and inverse of the magnetic field (A). 



Let us now turn our attention to the transverse Ising 
model (7 = 1) at finite temperatures. Classical corre- 
lation and QD for second nearest-neighbors are shown 
in Fig. 6. Note that QD still increases with temperature 
but this effect is feeble when compared with the same be- 
havior in XX model. We also observe that the classical 
correlation is typically greater than the quantum discord 
for any temperature. This is the opposite scenario in 
comparison with the XX model. Remarkably, for a weak 
magnetic field, QD available in the XX model overcomes 
that of the transverse Ising model. So, applications of 
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FIG. 6: (a) Classical and (b) quantum correlations between 
second nearest-neighbors in the transverse Ising model as a 
function of temperature (kT) and inverse of the magnetic field 
(A). 



the XY model in QIS tends to offer more quantum corre- 
lation as a resource in the isotropic limit (the XX model) 
than in the Ising limit. 



an analytical expression for both classical and quantum 
correlations for spin pairs at any distance. Remarkably, 
we have shown that the quantum discord between far 
neighbors is able to characterize a QPT, even for dis- 
tances where pairwise entanglement is absent. This is a 
consequence of the longer range of nonclassical correla- 
tion (that is nonvanishing for a class of separable states) 
in comparison with the short-range behavior of pairwise 
entanglement. Concerning the thermal effect onto corre- 
lations, we have shown how QD can be increased with 
temperature as the transverse magnetic field is varied. 
Moreover, we have also outlined the dominance of the 
QD over classical correlation for the XX model in oppo- 
sition to the Ising limit. Generalization of these results 
for larger subsystems and further analysis of the long- 
range behavior of QD for quantum critical phenomena 
including spontaneous symmetry breaking (at zero tem- 
perature) are left for a future investigation. 



IV. CONCLUSION 

In summary, we have examined pairwise QD and its 
classical counterpart in the thermodynamic limit of the 
anisotropic XY spin- 1/2 chain in the presence of an ex- 
ternal transverse magnetic field. We have considered the 
system at both zero and finite temperatures, providing 
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